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A method of creep analysis of orthctropic circular 
cylindrical and spherical shells subjected to exisymmetric loads 
has been developed in the present -work. A general study of creep 
behaviour of cylindrical and spherical shells subjected to uniform 
internal pressure has been conducted for a wide range of the values of 
anisotropy coefficients and creep law exponent . Importance of 
thermal stresses has been illustrated by a particular example of 
cylindrical shell. Analysis includes determination of stress 
redistribution, strain rates, stationary state stresses. 

Application of reference stress technique has been 
extended to the analysis of shells. 

Details of finite difference technique as applied to the 
present problem are discussed and suggestions for efficient 
computer use are included. 
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CHAPTER I 


INTRODUCTION 

A large number of modem industries involve systems 
having cylindrical and spherical shells sub je cted to high tempera- 
ture and pressure. Nuclear, steam and gas turbine power plants, 
chemical industries, supersonic transport are a few of such cases 
where the design of cylindrical and spherical shells is of extreme 
importance. It has been found that creep phenomenon plays a domina- 
ting role in the process of defoimation of such shells. One of the 
primary aims of a designer of equipment subjected to creep, is to 
offset the possibilities of failure by excessive deformation. 
Dimensional changes with time may lead to serious problems and 
ultimately to failure of the system. 

Creep defoimation of shells, made of isotropic materials, 
has been investigated by a large number of researchers. The 
assumption that a material starts off, and remains, in an isotropic 
state cannot be always justified. Materials become anisotropic 
during working and heat treatment processes and continue to do so 
under load. Creep analysis of structural members made of ortho- 
tropic materials has been limited to thin and thick cylinders. In 
the present work creep behaviour of shells, made of orthotropic 
materials, has been investigated. 

Reference stress method for estimation of creep defoima- 
tion has been applied to simple structures like beams and cylinders 
by many investigators. In the piesent woik it has been attempl 
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to investigate the applicability of reference stress techniques 
for creep analysis of shells. 

1 . 1 EE VIEW OP PAST WORK 

1.1.1 Creep Under Constant Uniaxial Stress : 

A tensile specimen under constant stress will deform with 
time. This deformation depends on three main parameters, stress, 
time and temperature. The most general creep equation is therefore 

e c = f (<r, t, T) ( 1 . 1 ) 

A useful first approximation is to limit this general 
function to a law of the form 

<r c = f i ( ) f 2 (t) P(T) (1.2) 

The sep aration of functions f ^ ) and fg (t) is common 

in most of work on creep and appears to be generally acceptable for 

'V 

the purpose of calculations on components.. 

The function f ^ ( <T~ ) has been chosen in many different 

(l)* 

ways. Kennedy gives a full summary . The most common forms are 
given below . 

Norton ) = L tf"' 12 

Mcvetty = A Sinix . ( <J~ / &- q ) 

m 

Garof alo f Air) = A Sinh (<r* / g-J (l.3) 

m m 

Johnson = J) {f . . + 

where L, A, D^, Dg, m, m., m 2 are material constants. 

* lumber within brackets refers to reference number. 




The most commonly used function is the power law 
attributed to Horton. She reason for its popularity is its simpli- 
city in stress analysis. 

A large number of alternative expressions are suggested 
for description of function f (t)^\ In practice, the time .dependence 
of creep in a complex alloy can be established only by curve fitting 


of experimental data. In such an emperical approach complex actions 
associated with shape of f 2 (t) are lumped together into one simple 
explicit time function. The following time functions are of 
practical use. 


Andrade f 2 (t) = (l+ exp (k^t) - 1 


.1/3 


1 

l 


/ (1.4) 


f 2 (t) = F t' 


n 


n. 


Bailey 

G-rahsm and Walles f 2 (t) = %. t^ 1 
b, k^, F, n, a^, are constants. 

Temperature has two effects on creep strain rate. Pirstly, 

a change in temperature has imm ediate effect on material constants. 

Secondly, temperature encourages changes in material structure. 

Some of important relations are reproduced here from reference 1 . 

C- c = A exp (- Q/RT) ) 

e c = A [t exp (-Q/RT)] J' (l.5) 

£ = A T exp (-Q/RT) J 

c 


where Q is the activation eneigy, R is the Boltzman's constant and 


A is a constant 
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1.1.2 Creep Under Multiaxial Stress : 

The first step in creep analysis under multiaxial stress 
system is to determine the combined effect of a multiaxial stress 
system to produce a certain magnitude of strain. This has led to 
acceptance of twin quantities, equivalent stress and equivalent 
strain increment, defined,' in general as 


cr = or (cr , cr , 


e 

A£- = 


x 7 

A ( A, 


or ) 

Z 7 


( 1 . 6 ) 


These quantities are usually chosen to reduce to correct values 


of stress and strain increment in the uniaxial case. For example, 

m 

an m-power relationship is = L <r At. By analogy 

with plasticity laws effective stress cT q is commonly assumed to 
take the same form as von-Mises yield function. For isotropic 


materials this is given by 

1 


cr _ 
e 


2 


m 1/2 


( <r - o- ) 2 4 ( cr- - <r ) 2 +( cr - cr f\ ( 1 # 7) 
x y J 7 z' x W J ' 


The next step is to determine principal strain increments. 

Extending the analogy with plasticity, equations similar to Prandtl- 

( 2 ) 

Reuss flow rule can be written as given below . 

= f 1 ( <r e ) K (1-8) 

( 3 ) 

Bhatnagar and Gupta ' derived constitutive equations for 

creep behaviour of orthotropic material. The form postulated is 

similar to the one for isotropic material. Effective stress and 

between 

relation -be hav iour strain rate and stress are 
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= — h -<r ) 2 +g (<r -cr _) 2 + h (<r -<r ) 2 ] 

a e ,/B“ b x y y 3' v x as'-J 


e v'2 
d'f „ 


dt^ = 2^“ C^ +H) ‘U - * r j - « <] eto ' 


(1.9) 


(1.10) 


•where ]?, G, H depend on material anisotropy. 

g^e equation (l.9) is based on assumption that creep behaviour is 

same in tension and c expression. 

(4) 

Ehutn sgar and Gupta used these relations for creep 
anaLysi® of thick walled cylinders. 

1.1.3 Creep Analysis of Structures ; 

A method for the determination of stationary creep of 
pressure vessel closures with constant stress rates using membrane 

theory of shells and povrer creep lav was developed by Gemma, Rowe 

(5) 

and Spahl 

The state of deformation and stress in a circular cylin- 
drical shell subjected to uniform lateral pressure was analysed by 
Bieniuk end Freudenthal ^ . A velocity field of deformed shell 
waS assumed with sane parameters which were determined with the 


odd of condition that power 


? dk E ,d € „ 

Yl ( ZSC ft f' 

^x dt~ + ~ L ~dt — * P tyj dx becomes minimum. 


Some approximations of the type 


d€ 

dt 


e 


m 


= f (t ,t) were used, 

,( 7 ) 


Sankamarayansn ' presented a method for stationary s tat e 
creep analysis of circular cylindrical shell under axially 
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symmetric load. Hie analysis was based on "square" interaction 

criterion between generalized stresses, namely, circumferential 

force and axial bending moment. It was assumed that creep strain 
is 

rate /a product of power function of a generalized strfess and a 
function of time. 

A simplified theory of stationary state creep of shells 
has been described by Rabotnov v . Hie real shell was replaced by 
a two layer model. A constant stress was assumed in each cf these 
layers. Thickness of layers and distance between them was so 
selected that behaviours of the real shell and its two layered 
model were similar in a moment-free state of stress and also in 
case of cylindrical bending. 

The methods described above may be used for stationary 
state analysis only and hence they do not give an estimate of 
redistribution of stresses* 

(q •iq') 

Penny ’ ' developed a method for calculating creep defor- 

mation in spherical shell (with axisymmetric openings) and clamped 
circular cylindrical shell. In this method, which determines the 
redistribution of stre.sses also, starting point is tie elastic 
solution. Stress rates and strain rates which depend on instant- 
aneous value of stress can be calculated, forward integration in 
time gives stresses after a small interval. Such a process may 
be repeated till stress rates throughout the structure become low* 

A method of analysis of creep bending in circular plates 

(l l) 

with temperature gradients was presented by Cozzarelli . The 
analysis was based on a power creep law in which the creep 



parameter was expressed as a function of temperature. Results 
were obtained for plates subjected to radially symmetric loads and 
temperatures. 


1.1.4 Reference Stress s 

Deformation rates in some structures which, obey power 
law can be estimated from the results of a single test conducted 
at a particular value of stress called reference stress. 

Marriott and Deckle studied stress redistribution which 

occurs during creep of a beam in pure bending and a pressurized 

thick walled cylinder, when elastic and primary creep effects are 

included. They found that magnitude of effective, stress during 

progressive stages of creep remained constant at some point which 

was named skeletal point. Both the components studied belong to 

a special class for which strain distribution is independent of 

the type of stress-strain relation. In such structures reference 

stress coincides with stress at skeletal point. The analysis has 

(l2) 

been described by Marriott in his review paper . 

(l3) 

Mackenzie ' show/ed that reference stress can be establi- 
shed by comparing the m-powrer solution w/ith that for m = 1. 

Results were obtained for beams in bending, pressurised thin 
walled cylinder and sphere, thick walled cylinder and axisymme- 
trically loaded circular plates. It may be remarked lhat 
application of this method is limited to structures for which 
analytical expressions for stationary state deformation are 


available 
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(14) 

Sim ' extended Mackenzie's approach to include the 
results of numerical solutions. Problem of bean in pure bending 

and of spinning disc mounted on rigid boss were solved. 

(15) 

Sim described a general method of obtaining reference 
stress from the dimensionless solutions of creep behaviour of a 
structure. The method requires stationary state solutions for 
two different values of creep law exponent. Reference stresses 
for a spinning disc mounted on a rigid boss and pressurised thick 
walled cylinder were obtained. 


1 .2 oBJ ECi ryiss ahd scope op the wor k 

All the earlier investigations into creep of shells 
assume that material is isotropic. Shis assumption is not a rea- 
listic one. Materials become anisotropic during working and heat 
treatment. Also materials may become anisotropic vfcen subjected 
to loads. 

In the present work a method of analysis has been 
developed to find creep deformation of shells made of orthotropic 
material. Uniaxial creep law of type 

id) ^ 

has been used throughout analysis. The main reason for choosing 
a power type law is that it has provided fairly accurate results 
in earlier investigations. At the same time analysis is made 
much simpler. .Moreover many materials like carbon steel, stain- 
less steel, aluminium alloys, nickel alloys and zirconium alloys 

(16 17) 

have been found to obey such a type of relation' 1 . 
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A constitutive law of the type given by equation ( 1 . 10 ) 
has been used to describe stress-strain behaviour in a multiaxial 
stress field. 

Cold working does give rise to difference in creep 

rates in tension and compression just as it produces preferred 

orientation. Ihe former can be removed by mild annealing. Hence 

it is perfectly possible, both in principle and practice, to 

have a material with identical creep characteristics in tension 

and compression and thereby satisfy assumptions inherent in 

equation (l.io). Zirconium alloys, for example, have been reported 

(l7) 

to follow a creep law of this type . 

Due to lack of data, change of anisotropic behaviour 
during creep has been neglected in the analysis. 

It has been shown that reference stress technique, 

■which greatly reduces the number of creep tests required to 
establish creep behaviour of a structure, can be applied to 
shells also. 

For' the class of materials described earlier the follow- 
ing problems have been solved. 

(i) Determination of stress redistribution, strain 
rates, steady state stresses and reference stress for (a) circular 
cylindrical shell with either fixed or simply supported edges 
(to) clamped spherical shell subjected to uniform internal 
pressure. 

This problem has been solved for a wide range of the 
values of stress exponent and anisotropy coefficients. 
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(ii) Determination . of stress redistribution, strain 
rotes and steady state stresses in a clamped cylindrical shell 
subjected to a uniform, pressure and axisymmetric linearly varying 
temperature field. 

Dhis problem has been solved for a shell operating in 
fast neutron flux and mode of zirealoy-2. 

ill the results aro presented in form of graphs, ill 
■|iie parameters associated with analysis are non-dimensional. 
However, the results presented are confined to a particular 
value of dimensionless length for cylindrical shells and of 
subtended meridian angle and radius to thickness ratio in case 
of spherical shells. 

Computer programs have been written for numerical evalua- 
tion of stresses and strains at different times during creep. 



CHAPTER II 


IMTHMATICAL EOEMUMTIOK gl METHOD OP SOLUTION 

In lii is Chester equations for elastic and creep analysis 
of shell problems stated in Chapter I have been derived and 
methods of solution are developed. 

2.1 ASSUMPTIONS 

(i) The material is homogeneous and ortho tropic. The 
axes of orthotropy coincide with the direction of principal stresses 
(ii) Power creep law stated below is applicable to creep 
deformation in uniaxial stress field* 

= B(t) p(t) cr m (2.1) 

> (iii) Von-Mises yield criterion and Prandtl-Reuss flow rate 

(3 ) 

for ortho tropic material are applicable. Stress-strain rate 
relations of type given by equation (l.io) are used*-. They- are -based 
on assumptions that creep behaviour in tension is similar to that 
in compression and creep strains do not cause any change in volume. 
Also creep is assumed to occur in any non-zero stress field. 

(iv) Displacements are small. 

(v) Hormals to the middle surface are conserved as such 
during deformation. 

(vi) Shell is thin so that z/a can be neglected in compa- 
rison with unity, This, in conjunction with (v), fofras Kirchoff- 

i 

love assumptions. 
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xs 

(vii) Structure/im elastic state at time equal to zero. 
This assumption is valid, for many practical cases. 

(viii) Load is constant. 

2.2 ELASTICITY EQUATIONS FOR CIRCULAR CYlIIIDRIC/i SHELL SUBJEC TED 
TO UNIFORM INTERNAL PRESSORE 

Sign convention is shown in Figure .2.1. 

2.2.1 Stress-strain Relations : 

For an orthotropic material stress-strain relations 


are 


( 21 ) 


<r * E h k _ + E. , 

x 1 xel 9el 


xel + E 2 fe 9el 


( 2 . 2 ) 


The above equation may be written in non-dimensional 


form as 


C = E 1 ^ xel * E 2 ^ eel 
JLm •»> M. 


E 


B 2 ^ 9el + % 


(2.2a) 


xel 


e 


E 


2.2.2 Strain-displacement Relations i 

These relations ere based on Kirchoff 1 s assumptions 
2 

c. _ dji _ d w 

x dx " z , 2 

dx 


(is) 


(2.3; 


9 


w 

a 


2.2,3 Stress-resultants and Stress Couples : 

Assuming Love’s first hypothesis 
h/2 


f 






(b) SPHERICAL SHELL 


G.2-1 COORDINATE SV 






(vii) Structure/in elastic state at time equal to zero. 
This assumption is valid for many practical cases. 

(viii) Load is constant. 

2.2 ELASTICITY EQUATIONS FOB GIBCULAE CYLIFDBICiL SHELL SUBJEC TED 
TO UNIFORM INTERNAL PRESSURE 

Sign convention is shown in Figure .2,1. 

2,2.1 Stress-strain Relations : 

For an orthotropic material stress-strain relations 


are 


( 21 ) 


^x ~ E 1 * xel + E ^eel 


<7“ s E C + e C- 
9 * 2 , xel *2 eel 


( 2 . 2 ) 


The above equation may be written in non-dimensional 


foim as 


< _ E 1 ^ xel * E 2 ^ Qel 

-r ~ # 


E 


E 2 ^ eel + ^ 


(2 ,2 a) 


xel 


e 


E 


2.2.2 Strain-displacement Relations 


These relations are based on Kirchoff's assumptions 
du 


,2 

d w 


dx 


dx 


(18) 


(2.3; 


e 


JL 

a 


2*2»3 Stress-resultants and Stress Couples 
Assuming Love’s first hypothesis 



If the ends of shell are closed. 


K -•§£ 
x 2 


Substituting equations (2.2) and (2.3) in the above 


equation, 


du _ o 1 w 

dx “ E 1 oi E 1 a 


(2.4) 


Based on Love's hypothesis , 

h/2 h/2 

i 


h/2 




© J ~ © 
-h/2 


dz , M 


z dz 


’ M e'J 


C - " zdz 


-h/2 


-h/2 


Substituting equations (2.2) through (2.4), the above 


equations are written in the following non-dimensional foim. 

E. 


_ a i w i 

19 _ e i * ■ * *0 

„ _ . J ML A sl/2 

9 /l ** ' / 

1 3E 
m © E 1 x 


(2.5) 

( 2 . 6 ) 
(2.7) 


Substituting equation (2.5) in equation (2.3), reduced 


strains are written as 
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2.2*4 Equilibrium Equations ; 

Equilibrium equations as given in reference 18 are 


dE 


x 


dx 


dx a 


* 

dx ~ Sc 


= 0 


= P 


> 


( 2 . 10 ) 


J 


Second equilibrium equation may be written as 


n_ 


'2 

< 


is. 


( 2 . 11 ) 


© «£ dv^ 

Eirst equilibrium equation and equation (2.4) are similar 
to each other. 

Combining equations (2.5) and (2.6) with the last two 
equilibrium equations, 


^4-^= o 

- d 1 


d4 

d V 


( 2 . 12 ) 


+ q = 0 


2.3 CREEP EQUAHOITS EOR CIRCULAR OYLIHPRICAL SHELL SUBJECTED 

TO tIHIEOHM IKTEMAL PRESSURE 

2.3.1 Stress-creep Strain Rate Relation j 

Assuming that strain rates depend on state of stress in. 
the structure, 
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= 

(a er ' 

V XX X 

d © 

xc 

at 

m-1 

= <r 

e 

a *ec 

m-1 

= r 

e 

dt 


A xx» \ 


ee © 


rt x@ ^~x 


^9 ) 


1/2 


(A X X “"i + A xe V P(T) B(t) 


^ er 0 + A, a <rj p(t) B(t) 


Vv , , 

/ (2.13) 


x© 


are called coefficients of anisotropy. A is taken, equal to uri J - 
for the sake of simplicity. A represents the ratio of creep ra+ - 
of uniaxial specimen with axis parallel to circumferential direc- 
tion of the shell to the creep rate of uniaxial spe cimen with axis 

parallel to longitudinal axis of the shell. A ^ measures the 

x© 

poisson effect due to creep in longitudinal and circumferential 
directions. 

Dimensionless time measure T is defined as 

t 


m-1 


¥ = cr Q P(T) J b(p) dp 

0 

t 

* r m 

_ J L . ? .( ? .) f b(p) dp 

* 


(2. 


, * \ 


cr Q 


The above expression is seen to be equivalent to the 
ratio of creep strain for uniaxial test at stress <r~^ to elastic 
strain due to stress <f^. 

Real time values can be calculated from Z by relation 

1 

(n-tl) T 


1 


„ * m-1 

L 105 ^0 


n+1 


( 2 . 15 ) 
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Non-dimensional form of equation (2 . 13 ) is 

1/2 

a XT) 

X© x ©' 



x = 

U r 2 + 

f 

£* H 


e 

• 

XX X 

m-1 

‘©© 

6 

V 

x = 

^xx 

r 

4* A 

xc 

X 

x© 

V = 

• 

> , = 

1 m " 1 Ow 

e ' ©e 


4 A 

2 

©c 

© 

X© 


( 2 . 16 ) 

(2.17) 

(2.18) 


2.3*2 Stress-resultant Hate and Stress Couple Rate 


(c = £ + €r 

x xel xc 


f = a + c 

v ^ ©c ©si 


i 

j 


(7 . . 


Equations relating rates of change of stress, strain and 
displacement are similar to equations (2.2) and (2.3 ). Combining 
such relations with equation (2.19), 


C~~ = E. (— — z r. 

x 1 v dx ,2 


au a 2 i 5 _ *_ (g A +E f ) ) 

' ^ a 1 xc h ©c' ‘ Vi 


cr 


dx 

w , ^ ,-du 


(2.2o) 


, 2 . 

d w 


e 


E 2l + E X- fe-%2 

dx 


) - (E„ € Q + E £ ) 
' '•2 ©c 2 / v x c / 


Since structure is subjected to constant load 
h/2 

= / K az = 0 


-h/2 


Substitute equation (2.20) in the above relation to get 
h/2 

-h/2 1 1 


dx h 


(2.21) 
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Love's hypothesis gives 


h/2 

h/2 

h/2 

r * 

f ’ * 

f 

1 j cr e dz ? M x = 

j (T x z dz ; M @ = 

<T~ e z dz 

-h/2 

-h/2 

-h/2 


The above relations are non-dimensionalized to the 


following form by using equations (2.20) and (2.2l). 

1 

* w 1 _ A 


n_ = - E __ 
S a o- 0 


j T 2^ 


( 2 . 22 ' 


1 


-1 

i dV f E 1 >,1/2 , E 1 is xx +E i/ A xe 

" ~ fe) - V4 


m x “ " s d n V 3E*' v T E * 


U 


^ d^ (2.23) 


321 


_ _ i a y dvfi %i/2 , *2*e© + B * A xe f u ]( „ 4 n 

« 2 E 1 d ^ ^ - 1/4 / ' M (2 * 24) 


2.3.3 Stress Rate and Strain Rate : 

Strain rate and displacement rate relation similar 

equation (2.3) is combined with equation (2.22) to get 
1 



Similarly equation (2.2) in teims of rates is 
substituted in equation (2.19) to get 


to 


(2.25) 
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a 

• i + (e„ /BJX e, a + e., a~ 

-5- _ Ax V_9 1 XX 2/ x6 T | 

* ~ *“ L 1 
X E E 1 

• „ _ Vee ^xo TI 

Q — -“O * *“ If U 0 

' E E 


/ (2.26) 

J 


2.3.4 Equilibrium Equations : 

Differentiate equation (2 . 1 1 ) to get its rate form 
given below. 


dq 

— + n = 0 

n e 


(2.27) 


Equations (2.22) and (2.23) are combined with, equation 
(2.10) after differentiating it with respect to time to get the 
following equilibrium equations. 


, 2 . 


* • 

- y = p 


where 


p • 

2 + 9 


(2.28) 


G 


1 


- *4* J T 2 S 

1 -1 


1 


1 


* 3 ( E \i- E„A + E_ A , f 

~ ~ * ° ( } ~ A j 

l -1 


c 2 v 3B„ 


E 


(2.29) 


2.4 ELASTICITY AMD CREEP EQUATIONS EOR A SPHERICAL SHELL SUBJECTED 


TO UNIEOHM INTERNAL PRESSURE 

Sign convention is shown in figure 2.1. Hie following 
equations can be. derived in the same manner as described in 
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sections (2.2) and (2,3). Details of derivation, therefore, have 
not been included. 


Equilibrium equations axe 


E 


E, 


<L 2 + ^)S-2* 2 @ - _ 2k 2 


2 ^ ■)©+ 2Jc 2 S 


(jf - 


E 


1 


* ' E 1 


I Z = 0 (2-30) 




(L 2 + j~) S - 2k 2 (h) = 2k 2 E c 


(L 2 -—-){§)+ 2k 2 S = 2k 2 G 


(2.31. 


where 


F c = - £ Cot 0 . J (V 1 - V 2 ) d^ + 
-1 


1 


j 


E, 1/2 




3E 


-1 

E, 


cot 0 <| 1 .. (A + A ^ 




-1 


lx’ -rf Bl )U 1^ ' (2 * 32) 


i ^ A ee ^ + e^ ^ u 2^ j 


1 

r 


d0 


E 


^ A xx +A x© E 1 ^1^ 


-1 


Stress couples and stress resultants are related to 
rotation ( 5 ) and shear force S according to the following 
equations. 

1 1 dS 

brirt 5 cot0 


(2.33) 
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m © 


J_h fg. 

12 a / 


(cot 0 (h) 


S 


E, 


d(H) 

d0 


E i a® e 


ilj I 


Reduced strains and non dimensional stresses are given by 
E v 


A 


n. 


v 


© “© E 1 “0 


n^ - £■ - t cot $ fE 
a > 


1 ^ i 

'V E 1 0 “ E 1 8 ' 2 a ’ a 0 

A 9 + EyA, A 0 

2 5 E 2 

© E 


> (2.34) 


r= 


X ^ A 


© 


E 


E 


1 


J 


Equations for calculating rate of change of stress resultants and 
stress couples are 


n © = 


dS 


1 


y 


djZf 

2k 2 

• 

Q 

cot 0 

o 


1 

h V 


“e ' ‘ 12 a 


E 


1 E 2 


Aqq + E ^ A x© 


E^dQi 

e 2 aF 

i 


E 


. 1 h E 1 / ^ 

m 0 = -l2* -5 ( T? 


J 

-1 


E 


2A 


E, 


cot 0 @) - 


( 2 . 35 ) 


. E .A +E. . A 0 
1 1 xx ^ x© 

I * 

4 E 


lM a{ ! 


-1 



Strain rates end stress rates may be written as 

Xi =*e-lf !i 0 + ^ j \ ^ - * ~A aot * © 

( . ** 


d © 


0 E 1 . '-© + - J W * a 1 W 


X© 


A© + ^ u// ^2 E 2 ^99 + 

'i “ * 


x© 


U, 


E 


E 


r- = E 
<>.J0 1 


. e a* + A x© 


u 


E 


E 


1 


(2.36) 


2.3 B IAST IOITY jgTD CHEEP BQU 1 TIGMS FOE L CIRCULAR CYLIHIRICfL SHELL 
SUBJECTED SO UNIEOBI IITIEHMAL PRESSURE AUD £XL SSMETBIO TEMPERATURE 
DISTRIBUTION 

The following stress-elastic strain relations are 
applicable 


cr~ — E,. _ + j ^ 


TP £ 

1 ' w xel‘ V ~ 9el 


(0 S e i + c< 2 s *> * 


n = E 2 «• ©el + \ t-xel- ( *2 E 2 + 0< 1 U 5 1 


"I 

j 

L 

( 

\ 

J 


(2.37) 


Equilibrium equations are 




■where 


■p _ J. _iL_ C . . ! ^ 
' £ c 2 d i V u - 

w *" “ > f ^ 2 

'* -I 


* _ 1 ( E* \ 1/2 E 1 il X£ + E i/ A xQ d f 

~ " 2 E * avj 4 


u 


$ 4€ J 


(2.40, 


Taking o( equal to unity stress resultants and stress couples may be 
related to rotation Y*" and shear force q. Thus, 


n e - 2 


m 


_ la 

n 

1 I2f (l±-\^ 2 

2 d *l V 


E, 


m e m x E 


1 

Strains and stresses may be calculated by using equations 


X 


E r 

— C 
x ~ E„ 


E 


Is 


<x 


ru -f 


E 1 ^ °<2 




% E h . 
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e 


E. 


I 
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Ae + E i // E 2 X x - (i + 
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(2.41) 


(2.42) 


Stress resultant rates and stress couple rates are given by 
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- - 


i E i/ d*/ / E i s 1//2 i E 2 A ee + E w A 


(-hi) 
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2 E 1 av l 3E 
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, IH 1 if, 
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Strain rates and stress rates are 
E 1 

X x = ~ eTJ - % + 2 J 7 1 d ^f 


dY-? /SeV^ 2 i 

d-n J \Eh J 


"1 


1 

Xe= £ Q + i I v 2 d ^ 

-1 


t- 

X 1 


. Ax + E v/ E i X © E 1 A XX + E t» A x9 , 
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E E 1 


V” 2 


X© + E t // E 2 X z E 2 A ee + E ^ A xe 


E 


E 


*2 


/ (2.44) 




2.6 BOUIMRY GOITEIglOUS 

The equations derived in the foregoing sections are to 
be solved subject to the following Boundary Conditions. 

(i) Circular cylindrical shell fixed at an edge (^= constant) : 

The necessary conditions are 

>6= 0 1 

V- 0 

f 

X 9= 0 

0 { 

(ii) Circular cylindrical shell s imply supported at an edge (> l =constant) 
= 0 


t= o 

t ^ o 


(2.45) 


m 


A 6 - 0 


(= 0 


(2.46) 


m 

x 

>9 


0 


U 


0 


0 
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(iii) Spherical shell fixed at an edge (0 = 0 ) 



= -0 
= 0 



0 


fc) = 0 1 

!• i >, o 
Ae = 0 > 

Also by symmetry, at apex of ihe spherical shell 



@ = 0 J 


IJIffiHiOl ® SOLUTION 


(2.47) 


(2.48) 


Finite difference method has been applied by many resear- 

(q 5 

chers to solve elasticity and creep problems of shells of revolution' 1 ’ 
An advantage of the method is that no approximations, other than those 
inherent in the finite difference method, are involved. 

To begin with, finite difference meliiod will be applied to the 
problem of circular cylindrical shell fixed at both the edges. Shell 
length is divided into n equal parts. Using parabola method of 
finite difference first and second derivatives may be written 


for shear force q 
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dq n _ jjn+1 ~ V-1 

d*-^ 2 A 

d 2 g Vj ~ 2 _S * Vjj 

a ^ ) 2 


(2.49) 


Similar expressions may be obtained for other variables also. 

Equilibrium equations may, therefore, be written as 


?or i = 0, 1, 2, n (fig. 2,l) 


q i-l " 2q i + q i-H 


aV. = 


° 


?i_1 “ 2 /i + ^i+i + ^ ~ 


2 . 


(2.50) 


Vi " 24 i + %+1 fi - A p oi 

. . . 2 2 . 

Y : * - 2f. +T-^ +A o. = a g . 

'1-1 'l 1+1 “l ^ Cl 


T^o 


y 


Similarly expressions for boundary conditions (2.45) are 
obtained (see Appendix A). Incorporating boundary conditions in 
equation (2.50), 2n simultaneous linear algebraic equations are 
obtained for each of the cases of elastic deformation ( X . - 0) 

'ia'feiL . _ , 

and oro ap deformation ( Z ?/ 0). These equations may be represen- 


ted by 


[a] j xj - |b] 

[a] [ij - 




(2.51) 

(2.52) 


Elements of the matrix and vectors , ^B^and 

|c| are given in Appendix A. 
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Similar matrix equations for simply supported cylindrical 
shell , clamped spherical shell subjected to uniform pressure and 
clamped cylindrical sheli subjected to non uniform temperature 
fields may also be obtained (see /ppendix A), 

It may be observed that coefficient matrix Qi.1 occurs both 
in the instantaneous equations (2.51 ) and rate equations (2,52). 
Vector depends on the loading only whereas vector | cj depends 
on current stress state and the load. 

Solution of equations (2.5l) and (2.52) is the complete 
solution of the problem. Steps involved in the solution are 
listed below. 


(i) Solution of the initial problem defined by equation (2,5l)» 

(ii) Evaluation of stresses, strains, stress resultants and 
stress couples by using vector J . 

6 

(iii) Substitution of Ihese stresses in expressions for # and 


G to get vector 


w 


(iv) Solution of the rate problem defined by equation (2,52). 

(v) Evaluation of stress rates, strain rates and rates of stress 
resultants and stress couples by using vector ^X j . 

(vi) Evaluation of new values of stresses and strains after a' 

small interval AT using expressions of Ihe type 2 = 

£ + 2 . * a Z . 

x z x z 

(vii) Steps (iii) to (vi) are repeated till steady state is 


attained 
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2.8 CALCULATION OP BEFBRBIK!E STRESS 

Hie main, difficulty in creep analysis even unde.n steady 
loads is the lank of a precise expression for the material behaviour 
viiich may be used conveniently. Using curve fitting techniques the 
data from constant load tests can be approximated by some simple 
expression. 

If structures are subjected to time-dependent loads the 
data from constant load tensile tests are completely inadequate. - ^ du]^6YS 

siwy*. 

q/2ietho d has been found of defining the material creep behaviour 
under variable loads. 

Hie concept of reference stress as a means of obtaining 
creep deformation In structures subjected to arbitrary load changes 
has been found to be promising. Hie creep behaviour of a structure 
can directly be linked to strain output from a single strain test 
conducted at reference stress. Reference stress approach avoids 
the need for defining material behaviour. 

It is possible to choose the unit stress, <3^, such that 
the parameters defining the stationary state behaviour of a 
structure are approximately constant over a range of values of 
stress index, m . Such a value of unit stress is defined as 
reference stress, <j~ 0 » 

It is unlikely that unit stress, <r* Q , first chosen to 
nondimen sionalize the problem will be exactly equal to the reference 
stress, o~ 0 ' It is desirable, therefore, to use dimensionless 
solution associated with an arbitrary value of unit stress to 
obtain parameters associated with reference stress. 
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Let the creep parameter of interest be strain S at 

c 

some point in the shell. Let this strain be reduced with respect 
to £ Q and denoted by A c • 

d A d K 

c _ i 6 c dt 

d Z £ 0 dt * dl 

substituting value of (dt/d T. ) obtained from equation (2.14) 
in the above equation 


If 


d£s c 

Tx~ 

~0 


P(T) B(t) cr o m 


lie 

dt 




(2.53) 

(2.54) 


dA 

then the dimensionless rate of parameter, £ associated 

with unit stress, cr~ o , is related to the dimensionless rate of 

d A c — 

parameter— associated with reference stress, <j~ q , by the 


following relation. 


* ILL 

^ c = ( ) A c (2.55) 

For m = O , & „ = A „ 

7 c c 

Hence for a particular set of results associated with an 
arbitrary unit stress, it is possible to calculate the constant 
using equation (2.55) and hence to determine reference stress 
using equation (2.54). The minimum number of points required to 
calculate the reference stress is two. However, when dealing with 
complex structures, solutions for more than two stress exponent 


values are desirable 
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It may be noted that except in certain, cases there is 
no single reference stress yfaich is applicable for all deforma- 
tion. rates in a structure. In general a particular reference 
stress is applicable only to a particular deformation parameter. 

It Eiay be seen from the definition of dimensionless time 
(equation. 2. 14) that its numerical value depends upon unit stress 
cr Q . It is possible to find a particular value of so that 

dimensionless time required to reach stationary stress, ~Z , is 

ss 

independent of stress exponent m. Ehis value of unit stress is 
called reference stress, associated with stationary state tine and 


is denoted by (Z ^ (=CK _ (T ) . It may be calculated by using 
OSS s s 0 

the following relation given in reference 15. 

El— 1 


z = (<x ) z 

ss ss ss 


For m = 1 , C ss =v ss 


( 2 . 56 ) 


2.9 HPMEB1CAL MEBI033S 

Matrix equations ( 2 . 5 1 ) and (2.52) have been solved by 

(l9) 

Choleski' s unsymmetric method . 

Is mentioned already, stresses, strains, and other 
quantities at any instant, are evaluated by equations of the 
following type. 


2. 


-2 +AI £ 


x 


Z +AX 


X 


T. 


4 x = A X ^ 

X t+ A T X X 


(2.57) 


/ 




Suitable choice of &T is very important, figure 2.2 
shows redistribution of stress at some point in tire shell 

starting Iron its elastic value X at Z= 0 to stationary state 

» 

? = £ ss * At "time T- 0 ? 2 q and X Q axe precisely known so 

* 

that the intercept of line X = X + X 0 7 wiUi the abscissa 
defines the point 7= Z . It is not desirable to choose 
^ X 0 / X 0 (= T q ) because Hi is would lead to absurdly hasty 

redistribution. Abetter choice is £T = crj/f; wli ere f ^ 1. In 
general 


A z = - 1 

f 


X 


I 


, f 1 (2.5s) 

An advantage of scheme given by equation (2.58) is that 
&~L chosen in "this manner by using local values of X and X Can 

4 

used throughout the computations. As time inereaseSj rate 
decreases and incremental time /1Z increases. A more general 
formula will be 

X m_ X „ iA 


- — 2 r minimum o 
f 


f ( 


jx 

a 

x 


x 


’ m ’ 


etc 


*)| (2.59) 


X 1U X X X ^ X 
It may be 'remarked that such a scheme will require exce- 
ssive computational effort. A better approach is to use equation 
(2.58) with a conservative value of factor f. 

When using equations (2.58) or (2.59) care should be taken 
to avoid those points in the structure where stresses (or strain 
etc.) -are very low. Very low stress may lead to incremental time 
approximately equal to zero. 
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Suitable choice of ^7 is very important. Pigure 2.2 
shows redistribution of stress at some point in the shell 

starting from its elastic value X at 7= 0 to stationary state ' 

* 

at Z = t . At time T- 0, 2 and X are precisely known so 
that the intercept of line X = 2 " 0 + X 0 7 with the abscissa 

defines the point 7 = It is not desirable to choose 

X X / X 0 (= T q ) because this would lead to absurdly hasty 
redistribution. Abetter choice is £ C = cr^/f; where f ^ 1 . In 
general 

a i I y 

j ) - 4 ^— , f 1 (2.58) 

Z 

An advantage of scheme given by equation (2.58) is that 

M 

chosen in ihis manner by using local values of X and X can 

4 

used throughout the computations. As time increases, rate XI 
decreases and incremental time increases. A more general 
formula will be 

1 . X m X , 

= — minimum of I ( — 7 — , -4 , ■ -t- — , , ... etc . ) j ( 2 . 59 ) 

f 1 £ ffi x Ax ! '^x 1 

It may be remarked that such a scheme will require exce- 
ssive computational effort. A better approach is to use equation 
(2.58) with a conservative value of factor f. 

When using equations (2.58) or (2.59) care should be taken 
to avoid those points in the structure where stresses (or strain 
etc.) are very low. Very low stress may lead to incremental time 
approximately equal to zero. 
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In addition to a suitable choice of incremental time, 
AT, it is necessary to decide when to stop computations. Static- 
nary sta^ may be defined based on one of the following criteria. 

(i) Strain rate is approximately constant 

(ii) Strain rate and creep strain rate are approximately 
equal 

(iii) Stress rate is approximately zero. 

In order to minimize computations without losing results 
of practical interest, it is suggested that computations be 
stopped when stress rate at a point, where strain rate is maximum, 
approaches a low value, this is necessary because of the asymp- 
totic nature of stress variation with time as steady s tate is 
approached. It is suggested that when stress rate at such a point 
is less than 0.01, the computations msy be stopped. 



results and d iscussion 


Results were computed for clamped circular cylindrical 
shell, simply supported cylindrical shell and clamped spherical 
shell subjected to a uniform internal pressure. Values used for 
the creep law exponent were 3, 5, 7, 9. Spherical shell problem 
was solved for an isotropic material whereas cylindrical shell 
problem was solved for the following values of the coefficients 
of anisotropy. 

°« 6 ^ A ee - 1,2 5 - °.3 2 4 ^ £ - o .6 

This range of properties covers a wide range of materials. A 
problem of clamped circular cylindrical shell made of Zircaloy-2 
and subjected to a uniform pressure and axLsynmetric linear tempe- 
rature distribution was also solved. 

The problems of cylindrical shell subjected to uniform 
pressure were solved for non-dimensional length of 80. The 
problems of spherical shell were solved for subtended angle 0^ 
equal to 30 degrees and ratio of radius to thickness equal to 50. 
The non-dimensional length of cylindrical shell subjected to linear 
temperature field was taken equal to 40. 

Unless otherwise stated, for all the problems E^ and Eg 
are equal. Ratio of to E^ is 0.3. 

3.1 SHELLS SUBJECTED TO UHIEOEM PRESSURE 

The main features of behaviour of shells are shown in 
Eigures 3.1 through 3.8 for clamped cylindrical shell, Figures 
3.9 through 3.13 for simply supported cylindrical shell, 

Eigures 3.14 through 3.20 for spherical shells. In all the cases 
unit stress is ^ . The results have been shown for nodes 
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starting from the shell edge to its middle. It is observed that 

(i) Axial stress at clamped edge of cylindrical shell is 
diminished from its initial value at Z = 0 to a stationary value 
at Z - ss by the process of creep (Figure 3.l). The higher Hie 
value of m, the creep exponent, the higher is the strain rate and, 
therefore, the quicker is the stress redistribution process. 
Similar stress redistribution is observed in the case of simply 
supported cylindrical shell (Figure 3.10) and clamped spherical 
shell (Figure 3.16). 

(ii) She effect of clamping in introducing stress concentration 
is experienced farther from the edge as compared to elastic case 
(Figures 3.3 and 3.4). 

Although creep process is powerful in diminishing peak 
stresses, the stationary state values of peak stresses are rela- 
tively insenitive to the values of creep law exponent m as is the 
distribution along the shell (Figures 3.4, 3.9 and 3.15). It is 
interesting to note that some skeletal points are present in ■ 
shells. Location of skeletal points is given in Table 3.1. 

The points where stress is independent of creep law exponent 
have also been included in the same table. 

The effect of redistribution of stresses on rupture 
times are tempting to speculation. If one were to use some 
sort of cumulative damage lav -in conjunction with rupture data 
from steady load tests, rupture times would be very much reduced 
from the values obtained by using stationary state stress._ 



TABLE 3.1 


3b 


Location of skeletal points along shell meridian. 


Type of shell 


T ~ 

,L_ 


Material 


T Position of jf 

-I £k^ls.tal_point_| 


Comments 


1. Clamped cylindrical 

Isotropic 


a 

Fig. 3.3 



= 0.7 

b 

Fig. 3.4 


V 

Anisotropic 

1! 

O 

• 

V>J 

b 

Fig. 3.4 



>| = 0.7 

b 

Fig. 3.4 

2. Simply supported 

Isotropic 

n = i.5 

b 

Fig. 3.9 

cylindric al 

Anisotropic 1 * 

i 

^ = 1.5 

b 

Fig. 3.9 

3. Clamped spherical 

Isotropic 

0 = 7-5° 

a 

Fig. 3. 14 



O 

VO 

CM 

11 

b 

Fig. 3. 15 


a Stress at stationary state is same as elastic value 

b Stress charges with time but its stationary state 
value is independent of stress exponent 

0 A ee =°- 6 > \ 6 =-°- 3 


(iii) Distribution of stress across shell wall is shown in 
Figures 3*2, 3.11 and 3.17 for clamped cylindrical? simply 
supported cylindrical and clamped spherical cases respectively. 
As expected, the distribution becomes nonlinear, unlike the 
linear distribution in elastic case. Higher the value of stress 
exponent, greater is the deviation from elastic solution. Again 
it may be observed that there exist skeletal points at 

(a) ^ = - 0.46 for clamped cylindrical shell 

(b ) Of = + 0.05 for simply supported cylindrical shell 

(c) ^ = - 0.90 for clamped spherical shell. 
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(iv) With reference to Figures 3.5, 3.7, 3.13, 3,19 it will be 
seen that rate of increase of strain is maximum at time T= 0. 

As Hie snell creeps, the stresses decrease and, therefore, 
strain rates decrease until they attain a minimum value at 
stationary state. 

It mjy be remarked that strains estimated by relation 

^ + 6" 35 t 

el ss 

will lead to reasonably accurate results for 5. 

Distribution of. circumferential strain in clamped 
cylindrical shell and spherical shell at different times 

<1 is shorn in Figures 3.6 and 3.18 respectively. 

(v) Variation of stationary state strain rates along shell 

meridian is similar to stress distribution at stationary state 

is 

(Figures 3.8, 3.12 and 3.20), Strain rate^ maximum. at points 
where stress is maximum. Thus in clamped cylindrical and 
spherical shells axial strain has maximum value &}; the edge. 

In simply supported cylindrical shells circumferential strain 
rate is maximum at the middle of length. 

(vi) Real time values can be calculated from ~Z by using 

equation (2.15). Time required to reach steady state may vary 

from a few hours to a few hundred hours (Table 3.2). It may be 

observed that a vessel operating for two shifts in a day will 

never attain a stationary state. Thus stationary state analysis 

(5 6 7 8 ) 

as advanced by many authors' 1 * ’ 3 ' may lead to non-conservative 


results 
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TABLE 3.2 


Time Required to Attain Steady State (a) 


Material 

f Type of shell f 

.1 L 

Temp 

(°c) 

1. b 

Clomped cylin- 
drical 

800 


Simply supported 
cylindric al 

800 


Spherical 

800 

2. c 

Clamped cylin- 
drical 

650 


Simply supported 
cylindrical 

650 


Spherical 

650 


I ra $ 0 \ Z l ss 

I _ _ Iteshs 1L ff JfesL-1. 


3 

400 

5.58 

31 

3 

400 

1.35 

8 

3 

400 

3.45 

19 

5 

500' 

1.25 

75 

5 

500 

0.95 

58 

s 

500 

. 3.11 

190 



(a) All tiie properties have been taken from reference 16. The 
creep law used in this reference is 



d £ 


dt 


10 


cr 


(f 

07 


18 Cr 8 Hi 0.45 Si 
E = 1.4 x 10 6 kg/ cm 2 


25 Cr 20 li 


E 


1.5 Si 
2 


1.4 x 10 kg/cm 


n 

■) 

0.4 Mn 0.1 C stainless steel 
2 

, = 50 kg/cm 

IMn 0.1 C stainless steel 
5 £T c7 = 350 kg/cm 2 


3.2 EFFECT OE AHISOTF-OPY 

A Q ~ and A 0 are two independent constants used to 
describe the anisotropy of an orthotropic material. As 
explained in reference 17, if A Qe < A^ ^then A^-^ -0.5. 
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Similarly if /” A^ , then <( - 0.5. Keeping this in 
view the following four cases were studied for simply supported 
cylindrical shell end- clamped cylindrical shell. Stress 


exponent in 

each 

case was 

taken as 

3. 

(a) 

A 96 = 

0.6 

"x0 

= - 0.3 


00 

A &© = 

0.8 

A xQ 

= - 0.4 


(c) 

A ee = 

1.0 

i; xe 

= - 0.5 

(isotropic Material 

(a) 

A e© 

1.2 


= - C.6 



The results have been 'plotted in figures 3.21 through 
3.26. The following observations are made. 


(i) Strain rates at stationary state are minimum for case (a) 
and maximum for case (d) (figures 3.21 and 3.22). Heferring to 
equations (2.17 and 2.18) it is seen that decrease in values of 


A q q and A^q results in decrease of creep strain rate. Comparing 


any two of examples under study it will be seen that an increase 


in the value of A and decrease in the value of Ag^ reduce 
strain rate. Thus it may be concluded that value of A ao is 


more important. 


(ii) As has been already remarked, higher -the creep rate, 
smaller is the value of maximum stress at stationary state 
(figures 3.23 "through 3.26). 


(iii) Higher the creep rate, quicker is "the stress redistribu- 
tion. Therefore, time to reach steady s tate is minimum for 
case (d) and maximum for case (a), (figures 3.23 and 3.24). 
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(iv) Variation of circumferential strain rate in simply supported 
cylindrical shell is shown in figure 3.27. Calculation of strain 
by superposition ox elastic strain and stationary values does not 
involve significant error. 

Spherical shell problem has not been solved for anisotropic 
material because degree of anisotropy of a spherical shell may vary 
along meridian and the shell tends to be isotropic at its apex 
when the shell properties are axisymme trie . 

Since anisotropy coefficients depend on texture of 
material, it may be possible to choose a fabrication process which 
develops a texture that gives minimum creep strain rate. 

3.3 THEMIX STRESSES 

A clamped cylindrical shell subjected to a uniform 

internal pressure and an axisymmetric temperature distribution 

varying linearly along the length of the shell was analysed. The 

results are shown in Figures 3.28 and 3.29. shell is made of 

zircaloy-2 which follows a law given by equation (2 . 1 ) wain 

^8390 

P (t) = ,0 x 26 x lO -1 '' e 1 (Reference 17) 

6 m = f 

Stationary state axial creep rate at the hot end was 

-7 -6 , • 

5.07 x 10 per hour as compared to 2.71 x 10 per hour at the 

cold end when normalized flux was taken as 2.7 and unit stress 

was 7500 psi. Maximum circumferential creep strain rate which 

-7 

occured at the middle of the shell was 0.69 x: 10 per hour. 

Creep caused by external loads leads to redistribution 
of stresses, drop in stresses at discontinuities being made up 



by increase in stress at other points. Since the thermaL stresses 
ere internally in equilibrium, they decrease throughout the whole 
-rolune of structure end tend, gradually, to disappear completely. 

To illustrate this, clamped cylindrical Shell subjected to uniform 
pressure and uniform temperature was analysed for tempera-fore rise 
of 280° c and compared with analysis for temperature rise of 30° c. 

In both the cases stresses at stationary state were the same. 
(Figure 3.30). 

3.4 EBFEEESrCE STRESS 

Figures 3.31 and 3.32 show the dependence of stationary- 
state strain rate and time l sg on stress exponent m. The fore 
of numerical results depends on foe magnitude of unit stress, 

o 2 ^) • i! - s m ®y seen^ there are particular values of unit 

stress at which foe strain rates and time T gs gr e approximately 
constant. Such a value of unit stress 

is the reference stress for foe particular parameter. Following 
points are observed. 

(i) Each of foe parameters, that is, axial strain rate at a point, 

circumferential strain rate at a point and time Z. has a different 

ss 

value of reference stress associated with it. 

(ii) Computed results for strain rates behave almost according to 
equation (2.55), the deviation being within + 4 per cent(Figure 3.51 
Therefore, it is possible to calculate , foe parameter defining 
reference stress, accurately. The accuracy cm be attributed to 
two factors. Firstly, veiy small change in foe value ofCK is 
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magnified m times when olo+fi-no- « 

P ' Secondly, since the dimension- 

loss strain rate xs conver^pfi"* - -r^ ^ 

eigS11 0 xo a fixed value, its numerical 

evaluation is relatively 

J senoitive to stationary state criterion 

used. 

(iii) Djjuensionless ttaes t ^ ao not confolm to eiJuation {2 . 56 ) 
to Ibe seme degree of ehcuracy ao strata ratos do to equation 
(2.55). The numerical methods used make it difficult to determine 
the stationary state time accurately. The results depend on 
otu.tionaiy staoe criterion used and on some programming ^xaM&t&rS. 
Scatter was found to be within + 15 and -4 percent. 

(iv) It may be seen from Figure 3.31 that reference stress for 
calculating the stationary state strain rate at a point in shell 
depends on anisotropy coefficients though it is independent of 
stress exponent m. Thus use of reference stress analysis to 
predict minimum strain rate by conducting one uniaxial creep test 
at reference stress is possible only when anisotropy coefficients 
are known. This is not a serious limitation since anisotropy 
coefficients can be related to texture of material • 

(v) Minimum strain rate obtained by reference stress analysis 
may bo used to gat strain at any time C according to equation 

A =A e(+ i X 

Strain calculated in this way constitutes a lower 
bound. Bie measured strain will be always greater than the 
predicted value using this method because of additional woik 
done as stresses redistribute from elastic to stationary state 



lues. 2hb difference, however, may not be nueh if structure 
is uncter creep for a sufficiently long tine. 

In addition, however, since the effect of stress 
redistribution is included in the original dimensionless solutions 
it is easy to obtain an approximate upper bound to deformation of 
suell by using an upper estimate of the stress exponent (Figure 3.33) 

3.5 NOTES Oil COMPUTATIONS 

'Hie computations involve solution of matrix equations 

of the typo 

M H = W 

Unsymetrical method of Cheloski was used. Step lengths 
of 0.1 end 0 . 5 ° were found sufficient for cylindrical and 
spherical smells respectively. Hie elasticity solutions were 
compared with analytical solutions given in reference 18 , Maxi- 
mum error in stress values was -0.2 per cent. Stress distribution 
in a shell subjected to uniform pressure should be symmetrical 
about the middle of length. Maximum deviation from sy mm etry of 
stress distribution was negligible for elastic solution aid 
1 percent for stationary state solution. Also by using 
Richardson extrapolation technique with step size equal tc 
one half of the step size suggested earlier, a maximum error 
of - 0.5 per cent was estimated. 

Simpson formula with ( m + 2) sections in half 

* • 

thickness was used to evaluate integrals like F and G . 

c c 
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Errors estimated by Richardson technique were found 

negligible . 


^cr oriental tine was calculated according to equation 

&Z. = 1 |(_|l_£)| for damped shells 

5 . x 

£~C = 4 j( - j for simply supported shells 
2- o 1 


Many values of factor f were tried. It was found that 
solutions no longer changed with f )>, 10. 


iill the: problems were solved on IBM 7044 computer. One 
problem requires 3 tc 15 minutes of computer time depending on 
the value of creep law exponent. Higher the value of exponent, 
more is the tine required. 



CHAPTER IV 


OONOHJSIONS 

In tne present work a method of creep analysis of 
cylindrical and spherical shells has been developed. The method 
can be applied to determine stress redistribution and strain 
accumulation in shells made of ortho tropic materials. Important 
conclusions from this work are listed below. 

(i) Creep leads to redistribution of stresses. High s tresses 
near discontinuities are reduced while stresses at other points 
of the structure increase. Since the effect of stress concentra- 
tion in shells is significantly less in creep than in elastic case, 
it may be inferred that membrane theory gives more accurate results 
in the former case. 

(ii) Stress redistribution is quicker and strain rates are- 
higher for higher values of creep law exponent. 

(iii) Ihen calculating rupture life of shelly stress redistribu- 
tion should be taken into consideration. Estimates based on 
stationary state values may lead to non conservative results. 

(iv) Strain rates are maximum at the time of loading and ihey 
decrease to a minimum value at stationary state. 

(v) Time taken by a shell to reach stationary state may vary from 
a few hours to a few hundred hours. Shells vhich are not in 
operation for long durations may never attain stationary state. 
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(vi) A reduction in the values of anisotropy coefficients A 

©9 

and resuuta in a decrease in stationary state strain rates 
rCid id 1 - increase in time required to reach stationary state. 


(vii) a d*. crease in the value of anisotropy coefficient A 


ee 


and on increase in the value of coefficient A within some 
bounds mjy also reduce stationary state strain rates. 


It any be possible to select fabrication processes to 
develop n texture which gives minimum stationary state strain 
rates. 


(viii) Thermal stresses diminish throughout the structure. 
Stationary 6 tatt stress distribution is independent of the 
temperature at which the shell operates. 

(ix) Inference stresses corresponding to stationary state V 
strain rate end time to stationary state have been calculated 
for 2 ^ a 9* With knowledge of anisotropy coefficients one 
uniaxial creep test conducted at the reference stress gives 
stationary state deformation at a point of interest. The 
necessity of evaluating the stress exponent can thus be obviated. 

( x ) Strain calculated by reference stress method is the lower 
bound to the stationary s tate strain. Its value calculated by 
direct creep analysis (using sn upper estimate of creep law expo- 
nent) Constitutes an approximate ijpper bound to- tire stationary 
state strain. 
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READ 1 » \'A' f , »‘AT HIT.A * ATZEE - 

READ it , DBS y* FACTS 

READ l f • » ALFA 

READ 1 , TFIKS T ,ILAS7, IMT "R, I ST P 
READ 1, IGAP1,IGAP2,IGAPD,IGAP > 

RF AD 2, pm AX 
PRiINT 3 
PRINT 4,NDIV 
PRINT 5, DELTA 
PRINT 6, EDS 
PRINT ? j FACTS 

PRINT 8, UMAX-' ' . . : , V E . . ■ .■■■.■EE,/..: : E . V E- ., ; >E : " X 

PRINT Hi ALFA ■ . - EE... ,/ / : eE ' E ■ ■. ... ' 'E . 

PRINT 9 ' ; ■: ■ v : ■ 

PRINT 29,Eli'H2,ENU,IMD -XiAXXi ATHITAiATZEE 
PRINT 1, IF I.RST.,.I'LA.STi INTER, ISTE.P 
PRINT 1, IGAflli IGAP2t IGA P3f I GAP4 
PRINT 2i EMAX 

ELEMENTS OF *-1A TP IX A ARE CALCULATED 
DIAGONAL TERMS ‘,P. r STORED AS VECTORS A»8»CrfD*E 


I L AS T 

INTER 

ISTEP 

I GAP 1 
I GAP 2 
I GAP 3 

I G AP 4 



o o n 


I PRINT =: 

INDEX* I NO \« « 

NUMB=NDI V*2-+2 

NEW=NUMB-1 

DO 3 V t 1 = 1 ,'fJHW, ? 

UI)=1, 

R ( I ) =Di'LTA^«2 

cm=-2, 

D( I) *f>. 

F ( I ) =1 e 

J=I+1 
AC J)=U 
B< j> =r . 

C( J) =-2. 

D ( J ) --DELTA* *2 
EC J ) =1. 

1X1/ CONTINUE 
At 2) =B« 

Rt NEW) =r. 

C( 2) *1. 

C ( NEW) =1« 

C C NU MB ) =2© 

D ( 2) =3* ■ 

D (NUMB )=( «, 

EC 2)=% 

€ ( 3 ) = 2<» 

EC NEW) =C: « 

ECNUMB)=U* 

CALCULATION OF ELEMENTS OF MATRICES L AND U SUCH THAT L*U=A 

Wt 1) =lc 
VCD =0. 

MU.) =0. 

Ntl)4 
N C 2) =% 

ZMR.I=0 « 

ZNRi*C. 

UC X) =C U) 


o o o o 
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Z = fi ( 1) 

ZMR=Z/ZUR 
M( 2) =Z MR 
Z = A( 3 ) 

ZNR1 =Z /ZUR, 

N( 3) =ZNR l 
DO 2 r :: I =2, NUMB 
W( I) = F.{ I ) 

ZWR=W( I ) 

Z = D ( I) 

Z VR.= Z- ZMR1+2WR 

V( I ) =Z.VR 
Z = C( I) 

Z U R= Z- Z MR *ZWR-ZMR*ZVP 

U( I ) =Z UP. 

IF Cl .EG* NUMB) GO TO 2'" i ■ 
ZMR1 =Z MR 

Z = 81 I) 

ZMR.= (Z~ZNR1*ZVR)/ZUP 
M ( 1+1) =ZMR 
ZNR,= ZMR1 

I F ( I •EQ« NEW ) GO TO 20*; 

Z = At I ) 

ZNR.1 =Z / Z UR 
N{ 1 + 2) = ZNR.1 
200 CONTINUE 


PLASTIC SOLUTION 

ITER ^COUNTER TO DETERMINE NUMBER OF STEPS 


ITR.=1 
I S ! TO P= 1 
PRINT V}. 

DO T 3. I =2, NUMB 
RHS( I) =L. 

Ifl CONTINUE 

RHS('l)=2.*DSLTA*(2*-*.;MU/fil)/ALFA 

■MEHRAfl';, , ..... , , 

CALL E CJU AT (M P H R k t- A L F « ) 

: i ' ^ 
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’/.izrfc AMD Mb; HR A Aft. INTRODUCED TO MATCH SUBROUTINE R.-QJIR M ; ‘T 
f ' JRIV,G is -MUMBrP OF DIVISIONS OF THICKNESS 

IN STR.f.SSCS AND STRAINS AND THEIR RATES FIRST SUBSCRIPT DENOTE 
NODE NUMBER AND SECOND SUBSCRIPT DENOTES POSITION ALONG THICK?) 

NRIN G1 = ( JNDSX+3)*2 
XINCR=2« /FLOAT ( NR I NG? 5 
NRING=NP ING1+1 
XX=E ?./ F } 

YY=FMJ/E1 

ESI AR, = F2"*ENU**2/E 1 
ZZ=S QR.T ( 3o *3ST AR /El) 

PP=E 1/ hST AR 
GQ=E 2/ ESTAR 
RR*ENJ/E2 

SS=( E1*AXX+BNU*ATZEE ) /ESTAR 

TT=( F2 *ATH IT A+fTNU* AT ZEE ) /ESTAR 

XK BAR 1 = ( FNU* AT H J T A+E 1*4 T 7EF ) / ( E 1* AXX+ EMU* A TZ E£ ) 

XKBAR2 = ( BNU* AXX+E 2*ATZEE ) / ( E2*ATHI TA+ENU*ATZEE ) 

TODEL=2.*DEUTA 
AA“2«/ ALFA 
B0=Y Y/ AL FA 
CC=XX/ ALFA 
Ml-N DI V+2 

DO 102 1=2, Ml ■■■■'■■ • .-0 ■ 'OX": • . 

NPHKI )=AA-tSH;-AIU 1 + 1 ) - SHEAR ( ) /T OD r L 

LAMDTl I ) =■' PH I( 1 )-?•?. 

XZI=*1. 

DO 1 f-B J = ; ,NRING ' .. ■ ' X.XX 

LAMDiXt I, J) =GC-YY«MPHI(I )-XZJ* ZZ* (GAMMA ( I +1 ) -GAMMA { I -1 ) ) /TODHL 
S I GM X { I , J ) = (, LA MD XU, J) + YY* L AM NT ( I ) } * P P 
SIGNTf I, J)=tLAMDT( I. ) +R-R.#L AM DX ( I , J ) )*QQ 
XZI=XZ I+XINOR 
103 CONTINUE 
1- 2 CONTINUE 
M I N 3 R= 1 
GO TO 2* •XI 
4002 TIME =0® 

PRINT 25, TIME 
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1 CALL DOTSf MRINIG, IN DRY, 4 XX, MHiTA ,XKRAR'i) 

CREEP SOLUTION 

LAST S T r P (4 rj)IS INCLUDED I ] BOTH ELASTIC AND CR’HP SOLUTION 
THIS ST: P CRTERMIM'-'S M0D1FI TD RHS FOR r QUATI ON A*Q = C 

V’- 1 MIN0R=2 

00 2L1 J=1 » NR I NG 

X = SI GMX ( 2t J 5 
Y=SI GMT { 2» J ) 

XYZ ! J ) =VH (Y ,X , INDEX* ATHITA ,A ;X) 

?J 1 COMT I N Or 

CALL QUAD! X‘YZ,NRING» XI ZGR ) 

XIZER = XIZER’/2. 

BEET A.= DELT A**?. 

RHS! 1) = B r :FT.A*FCDOT ! 2 ) +2 e *DF LTA.*XI ZER 
RHS! 2) *BS(~TA*GCDOT (2 ) 

Ml=2*NDIV+i 

L = 3 

DO nr 2 I=3.,Wl» 2 

RHS! I)=FCDOT{L )*B£ETA 

J=I + 1 

RHS! J ) =GCDOT ( L ) *BBETA 
L =L+ 1 

2* 2 CONTINUE 
MER.A=2 

CALL ‘ CUAKM HR A» XI ZER) 

CALC UL AT 10 I 1 IF 7R S " RATES AND STRAIN RATES 


Ml =>! 01 V + 2 

do ?. : j=:hmr 

NPOOT! I )=( SHEAR! I-D-SHEAR ! I+i) i/TODEL 

DO 302' J=1 iNRING 
X = SIGMX( It J > 

Y=SI GMT ( I »J ) 

V2( J ) = V11( Y,*X, INDEX, ATH I TA , AX X ) 

VI! J 5 = VU!X s fY» IHDE'XrAXXtATHIT A ) 

UI(J ) = U1KX.‘Y» INDTXt AXX » ATH IT A* XKBAR1 ) 



aan n no 
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3712 CO NT INUF: 

CALL CUAD( V2,NR INC ,T4DQ ) 

CALL QUADC VI, MR I JG »XADD ) 

LTDOK I ) =NPDOT ( J ) +TACD/2a 
XZ I = - 1, c 

DO 303 K=1 , NRI NG 

LX DDK I, K) =~YY*'!PnOT( I > - XZ J* ZZ* (GAMMA ( I +1 ) “GAMMA { I". ) ) / ( 2* »D'-LlA 
l+XADD/2, 

SXDOT ( I, K 5 = {,LXDOT{ I, K) + YY*LTDCT( I ) ) *PP~SS#UL ( K) 

ST DOT ( I , K ) = CLTOOT ( I ) +RR*LXDCT ( I »K> ) *QQ-TT*U2 ( K ) 

XZI = XZ I + XIN.CR 
303 CONTINUE 
301 CONTINUE 

CALCULATION OF SUITABLE INCREMENT OF TIME 


DTIME=10. 

DO 4fl 1 = 2, Ml 
DO ¥.■-:> J = 1,NRING 

IF ( A BS ( S IGMX (I , J ) ) o LT* Oo 1 ) GO TO 400 
IF( A8S ( SXDOT ( I , J ) } ®LTeEMAX) GO TO 40C 
RATI 0= S IGMX tl, J)/SXDOT( I, J) 

DT IMF1 =A8S ( RAT 10) 

I F ( DTI ME lo GToDTI ME )G0 TO 40 G 
DT IM E=DT IME1 

400 CONTINUE 

401 CQNT I MUX 

DT IMf= DT Ifi'T/FACTR 

TIME=T I ME +01 1 ME X'.I'-X- 

PAINT 23 


PRINT 25, TIME 
PAINT 1, ITER 


CALCULATION OF STRESSES AMD STRAINS AFTER D“LT A T I MS 





S I GMT ( I, J)=SIGMT{ I , J } +DT IMF -*$TDOT ( I , J ) 

f>i 1 CONTINUE 


CHECK FOR TERMINATION 


I F ( A BS ( SXO'OT ( 2 , NRI NG ) ) 0 GT* U BS )GO TO 2 /1 

PRINT 26 

PRINT 32 ? I TOR 

I ST0P=2 

MINI =3 

GO TO 6M 

2001 IF(ITER„LE* IFIR$T)MINI=1 
I F ( I TER* GTo IF I R ST ) MINI = 2 
I F ( I TER* GT* HAST ) M IN I =3 
I F(,I TER* GTo ITM AX ) GO TO 6C-8 
ITER=I TER+1 
IPRINT=IPRINT+1 
IF Cl PR INToGBo I NTER) I PR I NT= f 
IFCM IN IoNE, 1„ AND. IPR INToNE* 1) GO TO 4001 

60 1 PRINT 21 

00 602 1=2, ISTEP 

PRINT22, I ,>S IGMX ( I *1) , SIGMXC I , NR I NG ) , SIGHT! I » 1 ) » S I GMT ( I HIRING) , 
1LAMDXC I, 1},UAMDXU,NRING),LAMDT(I) 

602 CONTINUE 

J=ISTFp+i , ■■■V-;-' ■ ■■/ ■ "" , ' . ■: 

DO 603 I = J,Ml, I GAP A 


PRINT 22, I, SIGMXi I ,1),SIGMX(I , NR I NG) , SI GMT ( I , I ) , SI GMT ( I , NR ING) , 
1L AMD :-Cf I , ! 5 , LAMOM ( I , NRING ) , L AMDT ( I 5 
603 CONTINUE '//.■/.'■ . : V-/ ; ; H 

GO TOU. v,6 4) ,,'HnR 

6/;.4 go 
6'; 5 ll=m 

PR I M T 4 A 

DO 60 6 I=IGAP2,Ml, IGAP3 

PRINT 45,1, CSXDOT £ I » J ) , J=1 , NR ING , NR I N61 > » ( STDOT £ I , J ) » J = l > NRI NG, 
1NRING1 ) > (LXOOTU* J >, J=1 , NR ING, NRIMG1 ) , LTDOTt I ) 

606 CONTINUE 

DO 6* 7 I=IGAP1 ,M1 , LL 
PRINT 47, It ' T •' fc,BTMr - Y 


. 
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PRiIMT -a, I, (.SIGHT { I, J) , J = l, NRING) 

5f " 7 CONTI MU* 

GO T U( i j , 6 ' 9 ) , I STOP 
60 8 PRINT ;>7 
6C 9 STOP 

FORMAT SPECIFICATIONS 


1 FORMAT (I . IF) 

2 FORMAT (7H13« 6) 

3 FORMAT (* II NPUT DATA#//) 

4 FORMAT (KX.MTOTAL NUMBER OF DIVISIONS ALONG AXI S* * 43X, •"•=#, 15) 

5 FORMAT ( 1 X»*NON DIMENSIONAL LENGTH OF AN E L r MF NT*, 43 X, *=*, Fi4 e 8 ) 

6 FORMAT? If X»*STRF5S RATE ALLOWED AT STEADY STATE*, 44X, * = * , FI 4* 3 ) 

T FORMAT QOX.MFACTOP* BY WHICH STRESS/STRESS RATH SHOULD BE D9VIDLD T 
10 GET INCREMENTAL TIME =#,F14.8) 

8 FORMAT (I! X.*MAXIMUM HUMBER CF STEPS ALLOWED TO STEADY STATE*, 34X, 
1*=*, 15 ) 

9 FORMAT (/* PROPERTIES OF MATERIAL*// 

1* MODULI! OF ELASTICITY*, 24X,#CREEP LAW . ANISOTROPY 

2C0EF FIC IENTS*/* AXIAL* »14X,*CIRCUMFRANTIAL P 01 S S ON 

3 INDEX AXI AL*t i4X,*CIRCUMFRANTIAL POISSON*/) 

10 FORMAT (TF14o6) 

11 FORMAT (10X, ^MULTIPLYING FACTOR FOR REFPFNCE STRESS* ,37X ,*=#, F 12,, 5 ) 


TANG STRAIN 
OUTER 


21 FORMAT ( * MOD' AXIAL STRESS TANGl 

1 AXIAL STRAIN TANG STRAIN AXIA! 

2* OUTC INNER: 'OUTER 

3 !’ Ft "R-’O 

22 FORM AT ( 1 ) 

23 FORMAT? SX, ;.2 ( Ail-) ) 

2 5 FORM AT ( M X , * Tin- = > , F I 3„ ■ ) 

26 FORMAT?/* TEE SOLUTION HAS CONVERGED*// ) 

27 FORMAT?/# THE SOLUTION DOES NOT CONVERGE*/ /) 
29 FORM AT? 3? 5X^14* 6) ,?X,I 3,3? 5X,E14*6) ) 

31 FORMAT ( /* EUASTIC SOLUTION*//) 

32 FORMAT?/* TOTAL NUMBER, OF ITERATIONS =*,I5/> 
41 FORMAT?* INNER /OUTER, NODE, RAT 10 =*,2 15, FU.6) 

• n ..«A.n,ki ir.l.. jt'.! O.A'Ht.O. 1 1 C! 1, .. .L \ 


TANGENTIAL STRESS 
AXIAL MOMNT TANGMOMNF */ 
INNER OUTER 


42 FORMAT?* RATIO ~*,.11F1X 0 ? 

43 FORM AT ( 5 F 14.« 6 ) 




!■ IM 
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44 FDR M AT { *• NOD” AXIAL STRESS R 4TF*» 10X,*C IC UMFERENT I AL STRESS R. AT 
XtV X,» AXIAL STRAIN RAT”*,? X, *>C I R CUM FURORE NT J AL STRAIN RAT ) 

43 FORMAT ( IX, 14, 2X,1* 2 « <5 ) 

47 FORM AT ( * AXIAL STF.-SS =*•-! r ,£ 12e5/9Fi.% 5/ 9F1 ” ) 

48 FORMAT ! * TAN ST'-'.! SS =>-, I ,3 - 1 % 5 /<}«=!*, 5/ 9"i 3« A ) 

END 

I R FTC EQUAT 

SUBROUTINE ECU AT ( M r HPA, X I Z’ J R ) 

THIS SUBROUT I N b OR TERM! NTS TH ?. VALUES' OF SHEAR AND GAMMA 
MEHR.A=AN INBX WHICH SHOWS WHETHER THE PROBLEM IS ELASTIC OR 
IT IS OF CREEP 

XI ZER= VALUE OF INTEGRAL I ZERO (FOR CREEP ONLY) 

DIMENSION N 1 1 7 M , M ( ?. 70 ) , U { ?. 71) , V ( i 7 C > ) , W { 1 7 0 ) 

DIMENSION SHEAR ( 8 r > ) , GAMMA (85) 

DIMENSION RHS! 170) 

REAL N , M , 

COMMON/ Al/N ,'M, U, V,.W 
C0MM0N/A2/R.HS, SHEAR , GAMMA 
COMM ON /A3/ NO IV, DELTA 
COMMON/ A4/S1 » E 2 » EMU 
NUMB=2*NDIV-*-?. 

CALL MTR IX ( RHS, NUMB) 


1C l 
102 
im 

104 


OUTPUT OBTAINED FROM SUBROUTINE IS STORED IN RHS ITSELF 
NEXT STEP IS TO DISTRIBUTE VALUES GF RHS INTO GAMMA AND SHE .K 

V-CTJRS 

M1=MQI V+l OR'',;/ 

U = 3 

do i' : 7 i=r, m: 

SHEAR! I ) = RBS ! J ) 

GAMMA! I ) =R.HS ( J + 1) 

J=U+2 

CONTINUE 

GO TO !lU2 f |43 ) , MEH.R A 

SHEAR! 1 JsSHEAR ( 3 )-*?.« #Dr.LTA* 1 2®»ENU/E1 ) /XI ZER 

GO TO 104 

SHEAR! 1 ) =SHB AR ( 3 ) ** 2c *-DFLTA*XI ZER . 

SHEAR! 2 ) =f 
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GAMMA! i)=RHS(2 ) 

GAMMA! 2) =*;■«, 

MN=M1+1 
GAMMA! MM )= « 

SHF Aft! MM } =i. « 

MN=MI\1+ I - : 

SHF Aft! MN ) = RF S ( N U M p >» l ) 

GAMMA! MM ) =RFS ( NUMB ) 

RETURN 

END 

IBFTC DOTS 

SUBROUTINE DOTS ( NR IMG, I ND5 X ,A>’X, ATHI T A . XK3 AR 1 ) 

THIS SUBROUTINE DETERMINES THE INTEGRALS GCDOT AND FCOOT 
FROM THE KNOWN VALUES OF STRESSES 


DIMENSION S I GM X ( 8 5 » 2 B ) , SIGMT! 85 ,23) 
DIMENSION FQ!85),GC{ 35) ,FC DOT! 85), GCDOT (85) 
DIMENSION UK2 M,V?.( 2?) 

COMMON /AD/N'D IV , DELTA 
COMMON /A4/f;l , E 2, ENU 
COMMON / A 5 / S I GM X , SIGMT 
COMM ON /A6/FC DOT, GCDOT 
COMMON/ A7/U1/AIO/V2 
COMMON/ A 12/ X IN CR 
COMM ON /PARA /AT ZEE 
TQOFL=2«*DFliTA 
FORD EL =4* *D8LT A 

M=ND IV +2 . • 

DO' 101 1=2, M I ,;K~ R-. ■ '■ 

xzi=-l :r' 




n o o o non 


F r . ( I ) = F 

CALI ; IJAGt U , , M I H G » G ) 

GC( I ) = G 

i* 3 CDMiiNiyr ■ ■ ■■ ■ .. : ; 

~ =' _ ■!!)•■--' / 

i, p.i if,/ ) ••• ? : . .*•>;+ mu^atz”! ) n. 

rcnon 1 = ( fc ( '*-fc !".)} /T in i 
Fcnnif m ) = e 

GC m T ( ) ~ ’ >, { 0 C () *G C (=,})/ 0 PL T A 

GCD 3 K M ) ■ 

Jvi| | 3 ^ w ) y' B . 

nn i l i = ,m; . ■ 

pern r( i ) = r=ot i + )-f:. rr-i) s/ fowl 
gcdttc i )=;-;FKf gc< i-; )°gch + 1 ) ) /toosl 

IB CONTINUE' Z : ' I- BB . : '■ ' ZZ ■, Z: ■ ; .BbZ."' '■ p Z 

R.ET'JBM 

END 

IB FTC VXl 

FUMC TI ' )'! V \ l ( X , Y t I MO " • X ».R , S ) 

C 

C 'VALUATTB! OF CRCEP STR.AIN RATES FOR GIVEN STRESSES X, V 

I FX= AX I AL STRESS IT GIVES AXIAL STR.AI N D AT E 

IF ){=Z IRCUflFR AMTIAL STRESS IT GIVES Cl RCUMFRAMTI AL STRAIN RATE 


COMM OM / P AP.A/AT L 

VI 3 = S I GO,. ( X »'Y t X, R , S 5 * ( R *X+ATZEE*Y) 

R.FTJR.N 

-nd ;V g ! > z G 

1 B F T C U 1 1 

FUNC TI I' 1 ’ U ( BB PBSXf.R * $■? KBAR.) 

THIS F B 1 B i I ' B IE S U - UZ 

IF X* AXIAL STRGSS IT GIVES UI 

IF X = CT R,C U ■■ '1 F '* A N T I A L 'TR. SS ; IT GIVES U2 


cohmon/pab/./at r. ,z 
RJEAU < BAP. 

t j " ■ = 7 Z- ■■ f ;bb : :z V f :gb-Y) 

B-TUBi B'z 


■ i llifig i i ! In# ■ J ft 



IP FTC SIGN 

Fij't: rii." * t v ,r t ,y f t in- s) 

c 

c orrvi ar j- i'i tr a ctiv- s ir.~ss for given 

r 

COMMOg/P GP. ’ / ' T 7 
SIGN = r CRT ( P 

SIGN 1 * G I CM T jf ' 

P.FTJ R.M 
END 

I D FTC QUAD 

SUER.njTl' !!■ QUAD fF'.J 'Si I, M , RE S LT ) 

C 

C EVALUATION 'IF I IT GRXLS n Y SIMPSON FORMULA 

C INPUT I* A S' T nr VALU-S OF FUNCTION 

C 

DIME! IS I ON FUNS M ( a ) 

COMM ON/ A 12 /XI NCR 

R.FSL T = FUNSN L3 ) +PUNSNIN) 

M-N-3 • G ■ : 

on i n i 

A, Sl f= M ■'■+•• c FUNSMU 5 

4=1 + 1 

R,“SLT=r„ FSLT + 2**.BUWSMUr G 
ItJ'l CO NT I N U r . V G: V VyRr ■ G ' . G ;■ ;G . 

RiE--SUT.*'R E SCT+fe *FUNSN'< M+2) G 
■ G RJETL'= ' ! - .7 'V- /3* 

G J 

EG JO 

ib ftc aa.G 

coc- .■ o' a "a ■ f - < ' n : 

C .a G 

C THIS SUBROUTINE FINDS SOLUTION OF AX = P 

C C HOLY-S K I UNSYN METRICAL METHOD 

C OUTPUT V'CTOF P IS SCORED IN INPUT VECTOR P 

c numb - =TO T AL NUMBER OF UNKNOWNS 

DJtG-rGS ION ML •' ) » UO T r “ J »-V U.?G* ) * WU 7 f ) 

D I ME IS JO I PL’ "a • ; 

R.E'AL GJ * M : ' • : ■ ■A'' 1 :... -:G G ■■■.' MG: . . ■ ^ 


STRESSES X t Y 


ITSELF 



Be 1C 


COMMO 1/' ' i,. j, V, ;; 

,: ;OLJT ■ 0' i :>F L' • Y = P 

P( :>) = :; ( ' }• ) -:•- !( •':) 

DO 1 I ="■ 5 '■ j ! } '-1 3 

P( ] ) = P< i )*• M I 5 - P ( I ~ ■ >«■ i(I ) 

I i, C 1 HT I )!!' ■ , ■ 

NM = vlLM -o, 

P f 0 J *B ) = P( H j !■ 3 ) / U ( iiJt'BI 
P ( MM )= ( P ( ’M } 3( ) ,-.y( ) ) /UdlM ! 

no ? ’ i =p ,'j'i ■ 

J = NM ~ I +" 

P{ J) =( P( J ) R { J + '" } * 7{ J + ;. )-P ( J+ ,.)*W ( J + P. ) ) /UE J) 

?f t . CONT I 'jU* . . 

RETURN ' ■ ’ ■ 

END 

ENTRY 

\ 

' NOTE 

OTHER PROGRAMS USED IN THIS WORK ARE 

i l ) SR,? i:P ; ANALYSTS OF SIMPLY SUPPORTED CYLINDRICAL SHELL 
{2)CRd'P ANALYSIS OF CLAMPED SPHERICAL SHELL 


C.S: . MALY IF CL ‘.-IP C ,YU.DPIEfL SH L-L SU3J-CT i TO 

UNIFORM PRESSURE' AND AX! SYMMETRIC TEMPERATURE FIELD 
I M THESE PROGRAMS THOUGH THE EQUATIONS ARE DIFFRENT 
THE Y 1 T HOD OF COMPUTATIONS IS SIMILAR TO THE ONE FOR 

CLAMPED CYLINDRICAL SHELL 



